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Abstract 

The baryon asymmetry of the universe can be explained by the out-of-equilibrium decays 
of heavy right-handed neutrinos. We analyse this mechanism in the framework of a 
supersymmetric extension of the Standard Model and show that lepton number violating 
scatterings are indispensable for baryogenesis, even though they may wash-out a generated 
asymmetry. By assuming a similar pattern of mixings and masses for neutrinos and up- 
type quarks, as suggested by SO(10) unification, we can generate the observed baryon 
asymmetry without any fine tuning, if (B — L) is broken at the unification scale A GUT ~ 
10 16 GeV and, if m u ~ 3- 10 -3 eV as preferred by the MSW solution to the solar neutrino 
deficit. 
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1 Introduction 



The observed baryon asymmetry of the universe 

Y B = ^ = (0.6 - 1) • 1(T 10 , (1) 

s 

cannot be explained within the Standard Model, i.e. one has to envisage extensions of the 
Standard Model. Grand unified theories (GUTs) are attractive for various reasons and there 
have been many attempts to generate Y# at the GUT scale [|TJ. However, these mechanisms 
seem to be incompatible with inflationary scenarios which require reheating temperatures well 
below the GUT scale, the influence of preheating || on the baryon asymmetry requiring further 
studies. 

During the evolution of the early universe, the electroweak phase transition is the last oppor- 
tunity to generate a baryon asymmetry without being in conflict with the strong experimental 
bounds on baryon number violation at low energies ||. However, the thermodynamics of this 
transition indicates that such scenarios are rather unlikely ||. 

Therefore, the baryon asymmetry has to be generated between the reheating scale and the 
electroweak scale, where baryon plus lepton number (B + L) violating anomalous processes 
are in thermal equilibrium [Q , thereby making a (B — L) violation necessary for baryogenesis. 
Hence, no asymmetry can be generated within GUT scenarios based on the gauge group SU(5), 
where (B — L) is a conserved quantity. 

Gauge groups containing SO(10) predict the existence of right-handed neutrinos. In such 
theories (B—L) is spontaneously broken, one consequence being that the right-handed neutrinos 
can acquire a large Majorana mass, thereby explaining the smallness of the light neutrino 
masses via the see-saw mechanism 0. Heavy right-handed Majorana neutrinos violate lepton 
number in their decays, thus implementing the required (B — L) breaking as lepton number 
violation. This leptogenesis mechanism was first suggested by Fukugita and Yanagida [0] and 
has subsequently been studied by several authors PHHi. 

If one assumes a similar pattern of mass ratios and mixings for leptons and quarks and, if 
m v ~ 3 • 10~ 3 eV as preferred by the MSW solution to the solar neutrino problem, leptogenesis 



implies that (B — L) is broken at the unification scale | 12fl . This suggests a grand unified 
theory based on the group SO (10), or one of its extensions, which is directly broken into the 
Standard Model gauge group at the unification scale ~ 10 16 GeV. However, for a successful 
gauge coupling unification, such a GUT scenario requires low-energy supersymmetry. 

Supersymmetric leptogenesis has already been considered in refs. [p70] , |T~3|| in the approxi- 
mation that there are no lepton number violating scatterings which can inhibit the generation 
of a lepton number. Another usually neglected problem of leptogenesis scenarios is the nec- 
essary production of the right-handed neutrinos after reheating. In the non-supersymmetric 
scenarios one has to assume additional interactions of the right-handed neutrinos for successful 



leptogenesis (TT|. 
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In this paper, we investigate supersymmetric leptogenesis within the framework of the mini- 
mal supersymmetric standard model (MSSM) to which we add right-handed Majorana neutri- 
nos, as suggested by SO(10) unification. In the next section we will discuss the neutrino decays 
and scattering processes that one has to take into account to be consistent. In section |3| we will 
develop the full network of Boltzmann equations necessary to get a reliable relation between the 
input parameters and the final baryon asymmetry. We will show that by neglecting the lepton 
number violating scatterings one largely overestimates the generated asymmetry and that in 
our scenario the Yukawa interactions are strong enough to produce a thermal population of 
right-handed neutrinos at high temperatures. Finally we will see in section |] that by assuming 
a similar pattern of masses and mixings for leptons and quarks one gets the required value for 
the baryon asymmetry without any fine tuning, provided (B — L) is broken at the GUT scale 
and the Dirac mass scale for the neutrinos is of order of the top-quark mass, as suggested by 
SO(10) unification. 

In the appendices |A] and [FJ we introduce our notations concerning superfields and the Boltz- 
mann equations, respectively. The reduced cross sections for the scattering processes discussed 
in section can be found in appendix 0, while appendix [D] summarizes some limiting cases in 
which the corresponding reaction densities can be calculated analytically. 

2 The model 

2.1 The superpotential 

In supersymmetric unification scenarios based on SO(10), the effective theory below the (B — L) 
breaking scale is the MSSM supplemented by right-handed Majorana neutrinos. Neglecting soft 
breaking terms, the masses and Yukawa couplings relevant for leptogenesis are given by the 
superpotential 

W = ^N C MN C + iiH x eH 2 + HitQ\ d D c + E x eL\iE c + H 2 eQ\ u U c + H 2 eL\ u N c , (2) 

where we have chosen a basis in which the Majorana mass matrix M and the Yukawa coupling 
matrices A^ and A; for the down-type quarks and the charged leptons are diagonal with real 
and positive eigenvalues. 

The vacuum expectation values of the neutral Higgs fields generate Dirac masses for the 
down-type quarks and the charged leptons 

vi = (Hx) ^ =$> m d = \ d vi and m x = A; v 1 , (3) 

and for the up-type quarks and the neutrinos 

v 2 = (H 2 ) ^ =>- m u = \ u v 2 and m D — X v v 2 . (4) 
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The Majorana masses M for the right-handed neutrinos, which have to be much larger than 
the Dirac masses m D , offer a natural explanation for the smallness of the light neutrino masses 
via the see-saw mechanism ||. 

To get a no n- vanishing lepton asymmetry, one needs non-degenerate Majorana masses Mj. 
Then the scale at which the asymmetry is generated is given by the mass Mi of the lightest 
right-handed neutrino. Hence, it is convenient to write all the masses and energies in units of 
M 1; 



a.; 



Mj\ 2 s , Mi 



■ / •— \ , r i ■ % — T7o and z = — — , (5) 

J \Mj ' M\ T w 

where Mj are the masses of the heavier right-handed neutrinos, s is the squared centre of mass 
energy of a scattering process and T is the temperature. 

2.2 The decay channels of the heavy neutrinos 

The right-handed Majorana neutrinos Nj can decay into a lepton and a Higgs boson or into a 
slepton and a higgsino, while their scalar partners Nj can decay into a lepton and a higgsino 
or into a slepton and a Higgs boson (cf. fig. [[]). The decay widths at tree level read ||T3f 

~r N . := r(Nj ;^I+h) =r(N j -^V + h) 

r \ / - +\ M4 {m) n m n )ii , , 

= T( Nj -+l + H 2 )= T( Nj -> I + Hi) = ^ 1 ° vl D) " , (6) 

:= r(Nj^l + H 2 )^T(Nf^l + h) 

= r(N^f + Bl)=r^^i + h) = ^^^. (7) 

According to eq. ( P-5|) , the reaction densities for these decays are then given by 



(2) _ Mf {m) D m L 



_ 9/v W _ —1 V"D- b D)33 a 3\/ a 3 v f r— x (R s 

7iV, - 2 7 ^ - ^ ^ — K!(V%) • (8) 

All these decay modes are CP violating, the dominant contribution to CP violation coming 
about through interference between the tree level and the one-loop diagrams shown in fig. [1]. 
The CP asymmetries in the different decay channels of Nj and iVJ can all be expressed by the 
same CP violation parameter Sj, 

T^Nj ->T+h) -T^Nj ->P + h) r( y N j ^l + H 2 ) -T(Nj + H ] 2 ) 



£ 



' ~ r ( Nj ^i + h)+r(N j ^l t + h) r(N d -> l + H 2 ) +v(Nj ->7 + H\ 

r(ivf -+i + h)- t(n] ^i + h) r(ivf -> r+ h 2 ) - r(ivf + h\ 



r(Nf ^i + h)+ r(N$ ->i + h) t(n? ^T+h 2 ) + r(Nf ->r+H. 
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^2 (mlm^jj njtj 



(9) 



with = \/x 



In 



1 + x 



.r 



x — 1 



for i>1 



.r 



Here n is the flavour index of the intermediate heavy (s) neutrino. This result agrees with the 



one in ref. [13] and is of the same order as the CP asymmetry in ref. [10 
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Figure 1: Decay modes of the right-handed Majorana neutrinos and their scalar partners. 
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Figure 2: Contributions of the scalar potential to the decay width and the interactions of a 
scalar neutrino. 



With Ej we can parametrize the reaction densities for the decays and inverse decays in the 
following way 



7 



7 (iV/ 



>y(N$->l + h) = 1 (l + H 2 





=\ 


[1 + e j ) lNj 




=\ 


[I - Ej)^ N . 


>Nj) 


=\ 


3 




=\ 





Additionally, the scalar potential contains quartic scalar couplings, which enable the decay 
of N? into three particles via the diagram shown in fig. |^a. The partial width for this decay is 
given by 



64vr 2 vl 



r^:=r(Nf^T+U^ + ^) = Z^^^ with „„ 



Tr(AtA u ) 



47T 



(10) 



and the corresponding reaction density reads 

(3) _ 3a u M 1 4 (m) D m D )jj a jy /a] 



T 

1 J\ 



1287T 4 



(11) 



Since the Yukawa coupling of the top quark and its scalar partner is large, a u can be of order one. 



But even then 7— is much smaller than 7—. Hence, the three particle decays give only a small 

3 3 
correction, which we have taken into account for completeness. However, we have neglected 

the CP asymmetry in this decay which comes about through similar one- loop diagrams as in 

fig. E- _ 

The dimensionless squared total decay widths of Nj and Nj are then finally given by 



,( 2 ) 



aj (m) D m D ) 



jj 



Mi J 



167T 2 
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Mi 
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2 r 



.1.1 



J 



16tt 2 



1 + 



167T . 



(12) 



(13) 
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Figure 3: L violating processes mediated by a virtual Majorana neutrino or its scalar partner. 



The vertex in fig. [|a also gives 2 — > 2 scattering processes involving one scalar neutrino, like 
iVJ + / — > q + £/ c (cf. fig. Qb). The reduced cross section for this process reads 



0-22, {x) = 3a u 



4 



X 



(14) 



For the processes + q* — > I + f/ c and iV? + f/ c ^ — > I ^ + q, the corresponding back reactions 
and the CP conjugated processes we find the same result. The corresponding reaction density 



can then be calculated according to eq. (|B.6| ). One finds 



16tH 



3a u 
An a,jZ' 



InAz) 



(15) 



Hence, 722 will be much larger than 7^. for small ajz 2 , i.e. for high temperatures T ^> Mj. 



Together with similar scatterings which we are going to discuss in section |2.4j , these processes 
will therefore be very effective in bringing the heavy (s)neutrinos into thermal equilibrium at 
high temperatures where decays and inverse decays are suppressed by a time dilatation factor. 




Figure 4: Diagrams contributing to the lepton number violating scatterings via heavy sneutrino 
exchange. 

2.3 Lepton number violating scatterings mediated by the right- 
handed neutrinos 

Using the tree level vertices from figs. [I] and as building blocks we can construct lepton 
number violating scatterings mediated by a heavy (s)neutrino. Although of higher order than 
the tree level decays, these diagrams have to be taken into consideration to avoid the generation 
of an asymmetry in thermal equilibrium (cf. ref. [14||). In this section we will only mention the 



different processes which have to be considered. The corresponding reduced cross sections can 
be found in appendix and the reaction densities are discussed in appendix [D|. 

By combining two of the decay vertices (cf. fig. [I| and fig. |2]a) one gets the processes that we 
have shown in fig. and the corresponding CP conjugated processes. We will use the following 
abbreviations for the reaction densities 



7 « = 1 (l + h <-> P + h) , 7 i 2) = l(l + H 2 ^l + H\ 

if = 7 (T+ 1 ^ 7 + Hi) , T i 4) = j(l + h P + Hi 

7 (5) = 1 (7 , U. ^ P L 7Tc l *\ .,( 6 ) 



(l + H 2 ^V + U c + q) , ^ = 1 (l + H 2 ^l + h 



In' = 7 



The contributions from on-shell (s)neutrinos contained in these reactions have already been 
taken into account as inverse decay followed by a decay. Hence, one has to subtract the 



contributions from real intermediate states to avoid a double counting of reactions |L4 . 

From the scattering vertex in fig. |2|b and the decay vertices we can construct the following 
processes 
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7 (Z 1 + g <-> / + + # 2 ) = 7 (Z * + f/ c <-> / + g f + # 2 
7(7* + g «-> I + £ + f/ ct ) = 7(7* + f7 c «-> I + £ + g 



In fig. |] we have shown one typical diagram for each of these reaction densities. Again, these 
diagrams have on-shell contributions which have to be subtracted, since they can be described 
as decay of a sneutrino which has been produced in a scattering process. 



In ■ 



In ■ 



h 

h 
I 



I 

N 

--<-- I 

-<—H 9 



N 

-<—H 2 



+ 



+ 



h 
I 



N 



N 



/ I 

s T 



In ■ 



I 



h 



N 

-<—H 2 



+ 



I "- J 



■Ho 



1N C 



7iv • 



if 2 -->--r-<-- T 
q --< 



# 2 -->--n s ^ / 



In ■ 



->--r-<--^2 



— > 



g 



+ 



/ v yr-<-~H 2 
A ✓ i 

~ Y \ ' --Tic 

- g 



In ■ 



I 



N 



->--if. 



7jv • 



I -->■ 
ff 2 — «■ 



/ 



In ■ 



~* 1 *" 



h 



Figure 5: L violating processes mediated by a virtual Majorana neutrino in the t-channel. 

Up to now we have only considered processes with a neutrino or its scalar partner in the 
s-channel. In fig. |5| we have shown a selection of diagrams without on-shell contributions. The 
corresponding reaction densities will be denoted by 



In 

IN 



-f(h + h^ V + l ] , . 



7 S s )= 7 (i + /^fi} + Hj) , 



7 



7l 16) = 7 



= T (ff 2 + gt f +V + U c )=i(H 2 + fF +V + q) , 

7i 17) = -y(l+l++h + H 2 ) , ^=^(1+^^1 + 1), 

7l 19) = + P ^h + q^ + U c1 ) =l(l + q^l + U c] + h) 

= -f(l + U c ^T+q ] + h) = <y(P + h^J + q ] + tT ct ) 

= j(q + h^l+T+U c ) =j(U c + h<->l + l + q 1 ) . 

At first sight one may think that these diagrams could be neglected, since they are suppressed at 
intermediate temperatures, i.e. intermediate energies x ~ a,j. However, they give an important 
contribution to the effective lepton number violating interactions at low energies and therefore 
have to be taken into consideration. 

2.4 Interactions with a top or a stop 

The Yukawa coupling of the top quark is large. Thus we have to consider the lepton number 
violating interactions of a right-handed neutrino with a top quark or its scalar partner. 




Figure 6: Neutrino- (s)top scattering. 
For the neutrino we have to take into account the following processes (cf. fig. |]) 

Itf = l(N j +T~q + U*)=i(N j +T++q + v) , 
7§ } = l(Nj + q^P + U°) = 1 (N j + u^P + q) , 
llf = l(N J + uJ ~f f + q )= 1 (N J +q^V + u) , 
7t ( f = jfa + l^q + u) , 

4f = j(Nj + u I + g) = 1 ( y N j + q <-> 7 + u) . 

At this order of perturbation theory these processes are CP invariant. Hence, we have the 
same reaction densities for the CP conjugated processes. 
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Figure 7: Sneutrino-(s)top scattering. 



For the scalar neutrinos we have similarly (cf. fig. [?]) 
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The quartic scalar couplings of the sneutrinos give additional 2 — > 3, 3 — > 3 and 2 
which can be neglected since they are phase space suppressed. 
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Figure 8: Neutrino pair annihilation. 
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2.5 Neutrino pair creation and annihilation 

The Yukawa couplings of the right-handed neutrinos also allow lepton number conserving pro- 
cesses like the neutrino pair creation and annihilation. 

For the neutrinos we have the processes depicted in fig. || 

jfa + Nj^l + l) , 
j(N f + Nj ^h + h) . 



7 



(i) 



N C N C 
3 



JV?-->- 



7 



(3)_. 

N C N C ' 
3 



NP-->- 



h 



h 



7 



(2)_. 

N C N C ' 
* 3 



N c - 

3 



■>--T->- 
.<__! (_ 



" I 



X 

N Z 



7 



(4)_. 
N C N C ' 



iVf" ->"-!--> "-#2 



i Z + 

. < __L_<__iJ 2 



iVK. ,# 2 



Figure 9: Sneutrino pair annihilation. 



For the scalar neutrinos we have similar diagrams and additional contributions from the 
scalar potential (cf. fig. |9|). We have the following reaction densities 



3 « 3 



l~~ c =l{Nf + Nf ~h + h) , 7 g~ = 7 (iVf + iV| t ~tf 2 + tf 2 t 



' JVPiV c 



7 



(i) 



iVF-->- 



N; 



■~ I 



+ 



AT c -->- 



TV,. 



v# 2 

J — <- / 



7 



(2) 



jV c -->- 



N; 



+ 



ATP-- 



■> - ■ 



Zi 



Figure 10: Neutrino- sneutrino scattering. 
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Figure 11: Example of a Lf and L s violating MSSM process. 

Finally, there are neutrino-sneutrino scattering processes (cf. fig. pH|), 

& = 7(^ + ^^1 + 0, 7^ = 7 (A? + - * + H2) . 

J i 3 i 

Such diagrams also give neutrino-sneutrino transitions like Nf + I *->■ Nj + I. These processes 
transform neutrinos into sneutrinos and leptons into sleptons, i.e. they tend to balance out the 
number densities of the fermions and their supersymmetric partners, but they cannot wash 
out any generated asymmetry. As we will see in the next chapter, the number densities of 
the neutrinos and the scalar neutrinos are already equal without taking into account these 
interactions, while the equality of the number densities of leptons and sleptons is ensured 
by MSSM-processes (cf. section |2.6| ). Finally, the dominant contributions to these neutrino- 
sneutrino transitions come from inverse decays, decays and scatterings off a (s)top which we 
have already considered. Hence, we can neglect these additional processes. 



2.6 MSSM processes 

In the MSSM the fermionic lepton number Lf and the lepton number stored in the scalar 
leptons L s are not separately conserved. There are processes transforming leptons into scalar 
leptons and vice versa. As an example we have considered the process e + e «-> e + e (cf. fig. III]). 
For large temperatures, i.e. s 3> m^, the reduced cross section for this process is given by [15 
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This translates into the following reaction density 
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where we have introduced the dimensionless squared photino mass 



(16) 



(17) 



(18) 



These processes are in thermal equilibrium if the reaction rates are larger than the Hubble 
parameter H . This condition gives a very weak upper bound on the photino mass, 



m ^2.5xlO»Gev( T5 J^ 
In the calculations we assume m~ = 100 GeV. 
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412 V10 10 GeV 



(19) 
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3 Results 



3.1 The Boltzmann equations 

Now that we have identified all the relevant processes we can write down the network of 
Boltzmann equations which governs the time evolution of the neutrino and sneutrino num- 
ber densities and of the lepton asymmetry. For the scalar neutrinos and their antiparticles it 
is convenient to use the sum and the difference of the particle numbers per comoving volume 
element as independent variables, 



Y 



N i N c 
J 3 



(20) 



Furthermore, we have to discern the lepton asymmetry stored in the Standard Model particles 
Y L and the asymmetry Y La in the scalar leptons. 



For the neutrinos Nj one has 
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Correspondingly the Boltzmann equations for the scalar neutrinos read 



dY 



3+ 



dz 




- 2 



7 g } + 7§ + 3 722? + 2 7 f + 2 7 f + 2 7 £> + ^ + 2 7 g» 



+ - 



IY Y, 



2 Y~ Y~ q 

N c I 

3 



Y< Y L 



y^cq ycq Itj 



f (5) 

7tV 



+ E 



Y Y 

1 3+ 1 t+ 

Y~ Y~ 

JV9 N c 
3 * 



Y Y 

Y~ Y~ 

JV<? jvc 
3 * 



(fc) 



N C N C 

k=l 1 3 



+ 



Y Y 

I 3+ 1 



Y~ Y eq 




(22) 



dY 



Y 



dz sH(M 1 ) y: 



( € + T§ + 37-, + 27? } + 2^ + 2 7 f) + 7 « + 2 7 Jf 

3 3 



+ 



Y 



Y~ q 
i 



+ 



1 i 



(3) _ 1 (2) _ p (9) _ 1 ^3+ (8) L 1 

7^s j 2 ICS^ I 2 

J J Arc \ JV C 



722, 



(23) 



14 




15 



The numerical factors in front of the reaction densities arise due to the change in quantum 
numbers in the corresponding scattering, e.g. processes transforming leptons into sleptons ap- 
pear with a relative minus sign in the Boltzmann equations for Y Lf and Y Ls . Furthermore, any 
reaction density is multiplied by the number of different processes (cf. section |2|) contributing 
independently to the Boltzmann equations. 

This set of Boltzmann equations is valid for the most general case with arbitrary masses of 
the right-handed neutrinos. However, if the heavy neutrinos are mass degenerate, it is always 
possible to find a basis where the mass matrix M and the Yukawa matrix X u are diagonal, 
i.e. no asymmetry is generated. Therefore, one has to assume a mass hierarchy for the right- 
handed neutrinos, which in turn implies that the lepton number violating processes induced 
by the lightest right-handed neutrino are in thermal equilibrium as long as the temperature is 
higher than the mass of this neutrino. Hence, the lepton asymmetries generated in the decays 
of the heavier right-handed neutrinos are washed out and the asymmetry that we observe today 
must have been generated by the lightest right-handed neutrino. We will assume that the first 
generation neutrino N\ is the lightest. 

Hence, in the following we will always neglect the heavier right-handed neutrinos as free 
particles. However, they have to be taken into account as intermediate states since they give a 
substantial contribution to the effective lepton number violating processes at low energies. 

The fermionic part Y L of the generated lepton asymmetry will be transformed into a (B — L) 
asymmetry by the action of sphalerons. But since MSSM processes like the one in section [2.6| 
enforce the relation 

Y Lf = Y Ls , (30) 
the total lepton asymmetry Yi = Y L + Y Ls will be proportional to the baryon asymmetry ||16|| , 



YB --{22N f + l3N H ) YL > (31) 

where Nf is the number of quark-lepton families and Nh is the number of Higgs doublets. In 
our model with Nf = 3 and N H = 2 we have 

Ys = -±Y L . (32) 

From eqs. ([!]) and ([3(]) we can infer the asymmetries that we have to generate, 

Y Lf = Y Ls = -(0.9-1.4)-10- 10 . (33) 

The additional anomalous global symmetries in supersymmetric theories at high temperatures 
have no influence on these considerations, since they are broken well before the electroweak 



phase transition |17 |. 
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Figure 12: Typical solutions of the Boltzmann equations. The dashed line represents the equilib- 
rium distribution for the neutrinos N\ and the solid lines show the solutions for the (s)neutrino 
number and the absolute values of the asymmetries in Lf and L s , while the dotted line repre- 
sents the absolute value of the scalar neutrino asymmetry Yi_ . The lines for and Y\ + and 
for the two asymmetries Y Lf and Y Ls cannot be distinguished, since they are lying one upon 
another. The hatched area shows the measured value tffi\). 

3.2 The generated lepton asymmetry 

Typical solutions of the Boltzmann equations are shown in fig. [12|, where we have assumed a 
neutrino mass Mi = 10 10 GeV and a mass hierarchy of the form 

a 2 = 10 3 , {w) D m D )22 = a 2 (m j D m D ^ u , (34) 

a 3 = 10 6 , {rn ] D m D )^ = a 3 (m^m^^ . (35) 

Furthermore we have assumed a CP asymmetry E\ = — 10~ 6 . The only difference between both 
figures lies in the choice of (rn) D m D 



n 

4 



Km D ) n 10- 4 eV for fig. 12a, 
mi := V „/ = < (36) 
M i [ 10- 2 eV for fig. ||b. 

Finally, as starting condition we have assumed that all the number densities vanish at high 
temperatures T 3> Mi, including the neutrino numbers Y Nl and Yi + . As one can see, the 
Yukawa interactions are strong enough to create a substantial number of neutrinos and scalar 
neutrinos in fig. Oa, even if Y Nl and Y"i + do not reach their equilibrium values as long as z < 1. 



However, the generated asymmetries 

Y Lf = Y Ls = -4 ■ 10- 10 (37) 
17 



are of the requested magnitude. On the other hand, in fig. |12| b the Yukawa interactions are 
much stronger, i.e. the neutrinos are driven into equilibrium rapidly at high temperatures. 
However, the large Yukawa couplings also increase the reaction rates for the lepton number 
violating processes which can wash out a generated asymmetry, i.e. the final asymmetries are 
much smaller than in the previous case, 

Y Lf = Y Ls = -6 • 1(T 12 . (38) 

In both cases a small scalar neutrino asymmetry Yi_ is temporarily generated. However, Yi_ 
is very small and has virtually no influence on the generated lepton asymmetries. 

Usually it is assumed that one has a thermal population of right-handed neutrinos at high 
temperatures which decay at very low temperatures T <C Mi where one can neglect lepton 
number violating scatterings. Then the generated lepton asymmetry is proportional to the CP 
asymmetry and the number of decaying neutrinos and sneutrinos [|TJ, 

Y L « ex [ Y£(T > M x ) + Yg(T > Mi) ] w ^ . (39) 

With ex = — this gives 

Y Lf = Y Ls « -2 • 1(T 9 . (40) 

By comparison with eqs. ( |3"7D and ( |3"8"D one sees that by assuming a thermal population of 
heavy neutrinos at high temperatures and neglecting the lepton number violating scatterings 
one largely overestimates the generated lepton asymmetries. 

A characteristic feature of the non-supersymmetric version of this baryogenesis mechanism 
is that the generated asymmetry does not depend on the neutrino mass Mi and (m\ } m D )ii 
separately but only on the ratio mi [ITIJ . To check if this is also the case in the supersymmetric 



scenario we have varied mi while keeping all the other parameters fixed. In fig. [13] we have 
plotted the total lepton asymmetry Yl = Y L + Y Lg as a function of m"i for the right-handed 
neutrino masses Mi = 10 12 GeV, 10 10 GeV and 10 8 GeV, and the CP asymmetry e x = — 1CT 6 . 

The main difference between the supersymmetric and non-supersymmetric scenarios 
concerns the necessary production of the neutrinos at high temperatures. In the non- 
supersymmetric scenario the Yukawa interactions are too weak to account for this, i.e. additional 
interactions of the right-handed neutrinos have to be introduced. This is no longer the case 
here. The supersymmetric Yukawa interactions are much more important, and can produce a 
thermal population of right-handed neutrinos, i.e. the same vertices which are responsible for 
the generation of the asymmetry can also bring the neutrinos into thermal equilibrium at high 
temperatures. However, these lepton number violating processes will also erase a part of the 
generated asymmetry, hereby giving rise to the m% dependence of the generated asymmetry 
which we shall discuss in detail. 

First one sees that in the whole parameter range the generated asymmetry is much smaller 
than the naively expected value 4 ■ 10~ 9 . For low mi the reason being that the Yukawa in- 
teractions are too weak to bring the neutrinos into equilibrium at high temperatures, like in 
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Figure 13: Generated (B — L) asymmetry as a function of mi for Mi = 10 s GeV (dotted line), 
Mi = 10 10 GeV (solid line) and Mi = 10 12 GeV (dashed line). The shaded area shows the 
measured value for the asymmetry. 

fig. p^a. For high mi on the other hand, the lepton number violating scatterings wash out a 
large part of the generated asymmetry at temperatures T < M\, like in fig. |T2"|b. Hence, the 



requested asymmetry can only be generated if mj is larger than ~ 10~ 5 eV and smaller than 
~ 5 • 1CT 3 eV, depending on the heavy neutrino mass Mi. 

The asymmetry in fig. depends almost only on rhi for small mi ^ 10~ 4 eV, since in this 



region of parameter space the asymmetry depends mostly on the number of neutrinos generated 
at high temperatures, i.e. on the strength of the processes in which a right-handed neutrino can 
be generated or annihilated. The dominant reactions are decays, inverse decays and scatterings 
with a (s)top, which all give contributions proportional to mi to the Boltzmann equations at 
high temperatures, 

- - ~ -z (2) ~ -z ( 3 ) ~ 

1"i K mi ' tr/w \ 7~ oc mi , — — 7±J oc mi , 



sH(Mi) ' 1 ' sH(Mi) ' sH{Mx 



-z 



sH(Mi) 7221 K mi ' 7t(;)(T >> M0 K mi • (41) 

For large mi ^ 10 -4 eV on the other hand, the neutrinos reach thermal equilibrium at high 
temperatures, i.e. the generated asymmetry depends mostly on the influence of the lepton 
number violating scatterings at temperatures T^Mi. In contrast to the relations of eq. ([HJ ) 
the lepton number violating processes mediated by a heavy neutrino behave like 

7 f i (xM 1 ^mJ 1 i = A,...,D (42) 



sH(Mi 
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at low temperatures. Hence, one expects that the generated asymmetry becomes smaller for 
growing neutrino mass Mi and this is exactly what one observes in fig. ITS. 



The lepton number violating scatterings can also explain the small dependence of the asym- 
metry on the heavy neutrino mass Mi for mi ^ 1CT 4 eV. The inverse decay processes which 
take part in producing the neutrinos at high temperatures are CP violating, i.e. they generate 
a lepton asymmetry at high temperatures. Due to the interplay of inverse decay processes and 
lepton number violating 2—^2 scatterings this asymmetry has a different sign compared to the 
one generated in neutrino decays at low temperatures, i.e. the asymmetries will partially cancel 



each other, as one can see in the change of sign of the asymmetry in fig. |12|a. This cancellation 
can only be avoided if the asymmetry generated at high temperatures is washed out before the 
neutrinos decay. At high temperatures the lepton number violating scatterings behave like 



— z 



sH(M 1 



7f L oc Mi ajm] 



D . 



(43) 



Hence, the wash-out processes are more efficient for larger neutrino masses, i.e. the final asym- 
metry should grow with the neutrino mass Mi. The finally generated asymmetry is not affected 
by the stronger wash-out processes, since for small rh\ the neutrinos decay late, where one can 
neglect the lepton number violating scatterings. 

This change of sign in the asymmetry is not observed in fig. |i~2]b. Due to the larger m\ 
value the neutrinos are brought into equilibrium at much higher temperatures, where decays 
and inverse decays are suppressed by a time dilatation factor, i.e. the (s)neutrinos are produced 
in CP invariant scatterings off a (s)top. 



4 SO (10) unification and leptogenesis 



In ref. [[L2| , it was shown that there is no direct connection between CP violation and leptonic 
flavour mixing at high and low energies. Furthermore, the implications of non-supersymmetric 
leptogenesis on the scale of (B — L) breaking and the light neutrino masses have been studied 
by assuming a similar pattern of masses and mixings for the leptons and the quarks. Here we 
are going to repeat this analysis for the supersymmetric scenario. 



4.1 Neutrino masses and mixings 

If we choose a basis where the Majorana mass matrix M and the Dirac mass matrix mi for the 
charged leptons are diagonal with real and positive eigenvalues, 



mi 



m e 
m M 
V m T j 



\ 




M 



/ Mi \ 

M 2 
\ M 3 J 



(44) 
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the Dirac mass matrix of the neutrinos can be written in the form 



m. 



V 



( m x \ 

m 2 
V m 3 / 



Z7+ 



(45) 



where K and £7 are unitary matrices and the are real and positive. 

All the quantities relevant for baryogenesis depend only on the product m) D m Dl which is 
determined by the Dirac masses and the three mixing angles and six phases of U. Five of 
these phases can be factored out with the Gell-Mann matrices Aj, 



U 



e" ' e"'°" e"'""" U\ e""°" e 



(46) 



In analogy to the Kobayashi-Maskawa matrix for quarks the remaining matrix U\ depends on 
three mixing angles and one phase. In unified theories based on SO(10) it is natural to assume 
a similar pattern of masses and mixings for leptons and quarks. This suggests the Wolfenstein 
parametrization |TB[ as ansatz for U\, 



1 - 



A 



A\ 3 (p-ir]) \ 
-A 1 - f AX 2 

\ A\ 3 (l - p - irj) -AX 2 1 / 



(47) 



where A and \p + ir]\ are of order one, while the mixing parameter A is assumed to be small. 
For the Dirac masses SO(10) unification motivates a hierarchy like for the up-type quarks, 



mi = b\ 4 rri3 mi = c\ 2 m^ b,c = 0(1) 



(48) 

We have mentioned in section that we also need a hierarchy in the Majorana masses to 
get a lepton asymmetry. We choose a similar hierarchy as in eq. 



Mi = BX 4 M 3 M 2 = CX 2 M 3 B,C = 0(1) . (49) 

Later on we will vary the parameters B and C to investigate different hierarchies for the right- 
handed neutrinos. 



The light neutrino masses read |12 



m. 



\C + e 4ia B 

c 2 \C_ + e Aia B 
BC 



A m VT + 0(X { 



X 2 m u +OX 



(50) 

(51) 
(52) 



We will not discuss the masses of the light scalar neutrinos here, since they depend on unknown 
soft breaking terms. 
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In section |3.2| we have seen that the lepton asymmetry is largely determined by the mass 
parameter m 1; which is given by 



mi 



c 2 + A 2 \p + ir)\ 
B 



A 77V 



i.e. mi is of the same order as the mass. From eq. (|^) we get the CP asymmetry 
3 B A 2 - 2 



Si 



A 



4 ™ 3 Im 



]lT C 2 + A 2 \p + 77?| 2 v\ 



( p _ ir j) 2 e i2(a +v ^3/3)l + Q (\ 



In the next section we will always assume maximal phases, i.e. we will set 



B A 2 | 


P + 777 


2 


c 2 + A 2 


I/O + ir] 


2 



A 



4 m_ 



+ O (A 6 ) 



(53) 



(54) 



(55) 



Hence, the lepton asymmetries that we are going to calculate may be viewed as upper bounds 
on the attainable asymmetries. 

Like in the non-supersymmetric scenario a large value of the Yukawa-coupling m-&jvi will 
be preferred by this baryogenesis mechanism, since Si oc m|/w|. This holds irrespective of our 
ansatz for neutrino mixings. 



4.2 Numerical results 

The neutrino masses (0)-(0) can be used to constrain the free parameters of our ansatz. The 
strongest hint for a non-vanishing neutrino mass being the solar neutrino deficitQ we will fix 
the mass to the value preferred by the MSW solution |19|| , 

77v ~ 3 • 1(T 3 eV . (56) 



Hence the parameter mi, which is of the same order as 777.„ according to eq. (|53|) , will be in 
the interval allowed by fig. |13|. 

The most obvious parameter choice is to take all 0(1) parameters equal to one and to fix 
A to a similar value as the A parameter of the quark mixing matrix, 

A = B = C = b = c = \p + irj\ ~ 1 , A ~ 0.1 . (57) 

The mass in eqs. ( |51~D and (pED then fixes the v e and v T masses, 

~ 8 • 10" 6 eV , m Ur ~ 0.15 eV (58) 

and mi reads 

mi ~ 3 • 10~ 3 eV . (59) 
SO(10) unification suggests that the Dirac neutrino mass 7773 is equal to the top-quark mass, 

m 3 = mt ~ 174 GeV . (60) 

-"Tor a review and references, see pPf . 
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Figure 14: Generated asymmetry if one assumes a similar pattern of masses and mixings for 
the leptons and the quarks. In both figures we have A = 0.1 and = rn t (a) and m 3 = mi, (b). 



This leads to a large Majorana mass scale for the right-handed neutrinos, 

M 3 ~ 2 ■ 10 14 GeV => Mi ~ 2 ■ 10 10 GeV and M 2 ~ 2 ■ 10 12 GeV , (61) 

and eq. (j55|) gives the CP asymmetry E\ ~ —6 • 10~ 6 . Integration of the Boltzmann equations 
yields the (B — L) asymmetry (cf. fig. |H|a) 

Y B _ L ~ 1(T 9 , (62) 

which is of the correct order of magnitude. It is interesting to note that in the non- 
supersymmetric scenario one has Yb-l — 9 • 1CT 10 for the same choice of parameters. 

Our assumption (|6"0"D, m 3 ~ m t led to a large Majorana mass scale M 3 in eq. (0). To check 
the sensitivity of our result for the baryon asymmetry on this choice, we have envisaged a lower 
Dirac mass scale 

m 3 = m b ~ 4.5 GeV , (63) 

while keeping all other parameters in eq. (|57|) fixed. The assumed mass ( |56D then requires 
a much lower value for the Majorana mass scale, M 3 ~ 10 11 GeV and the CP asymmetry 



E\ ~ —4 -10 becomes very small. Consequently, the generated asymmetry (cf. fig. |l4|b) 

Y B _ L ~ 10- 12 , (64) 

is too small by more than two orders of magnitude. We can conclude that high values for both 
masses m 3 and M 3 are preferred. This suggests that (B — L) is already broken at the unification 
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Figure 15: Generated lepton asymmetry if one assumes a similar mass hierarchy for the right- 
handed neutrinos and the down-type quarks. 



scale A GUX ~ 10 16 GeV, without any intermediate scale of symmetry breaking, which is natural 
in SO (10) unification. Alternatively, a Majorana mass scale of the order of 10 12 to 10 14 GeV 
can also be generated radiatively if SO(10) is broken into SU(5) at some high scale between 
10 16 GeV and the Planck-scale, and SU(5) is subsequently broken into the MSSM gauge group 
at the usual GUT scale ~ 10 16 GeV 0. 

In eq. fl4D| ) we had assumed a mass hierarchy for the heavy Majorana neutrinos like for 
the up-type quarks. Alternatively, one can assume a weaker hierarchy, like for the down-type 
quarks by choosing 

B = 10 and C = 3. (65) 
Keeping the other parameters in eq. (|57D unchanged fixes the v e and v T masses, 

m Ve ~ 5 • 10" 6 eV , m Vr ~ 0.7 eV , (66) 

and the mass parameter mi, 

mi ~ 10~ 3 eV . (67) 
Choosing the Dirac mass scale (|60D we get a large Majorana mass scale 

M 3 ~ 4 • 10 13 GeV Mi ~ 4 • 10 10 GeV and M 2 ~ 10 12 GeV . (68) 

From eq. ( p5| ) one obtains the CP asymmetry e-y ~ —6 • 10~ 5 . The corresponding solutions of 



the Boltzmann equations are shown in fig. 15. The final (B — L) asymmetry 



Y B _ L ~ 2 • 10^ 8 , (69) 
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is much larger than requested, but this value can always be lowered by adjusting the unknown 
phases. Hence, the possibility to generate a lepton asymmetry does not depend on the special 
form of the mass hierarchy assumed for the right-handed neutrinos. 

In the non-supersymmetric scenario one finds for the same parameter choice 

Y B _ L ~ 2 ■ KT 8 . (70) 

Hence, when comparing the supersymmetric and the non-supersymmetric scenario, one sees 
that the larger CP asymmetry in the former and the additional contributions from the sneutrino 
decays are compensated by the wash-out processes which are stronger than in the latter. 

5 Conclusions 

We conclude that the cosmological baryon asymmetry can be generated in a supersymmetric 
extension of the Standard Model by the out-of-equilibrium decays of heavy right-handed Ma- 
jorana neutrinos and their scalar partners. Solving the Boltzmann equations we have shown 
that, in order to be consistent, one has to pay attention to two phenomena which can hamper 
the generation of a lepton asymmetry. 

First, one has to take into consideration lepton number violating scatterings. These pro- 
cesses, which are usually neglected, can wash out a large part of the asymmetry if the Yukawa 
couplings of the right-handed neutrinos become too large. 

On the other hand the neutrinos have to be brought into thermal equilibrium at high 
temperatures. We could show that for this purpose it is not necessary to assume additional 
interactions of the right-handed neutrinos in our theory, since the Yukawa interactions can be 
sufficiently strong to produce a thermal population of heavy neutrinos at high temperatures, 
while still being weak enough to prevent the final asymmetry from being washed out. 

The observed baryon asymmetry can be obtained without any fine tuning of parameters 
if one assumes a similar pattern of mixings and Dirac masses for the neutrinos and the up- 
type quarks. Then the generated asymmetry is related to the mass and fixing this mass 
to the value preferred by the MSW-solution to the solar neutrino problem leads to a baryon 
asymmetry of the requested order, provided (B — L) is broken at the unification scale, as 
suggested by supersymmetric SO(10) unification. 

In supersymmetric theories there are further possible sources of a (B — L) asymmetry, e.g. 
it may be possible to combine inflation with leptogenesis (cf. refs. |^2j|). In this connection, 
possible constraints on the neutrino masses and on the reheating temperature from lepton 
number violating processes at low temperatures require further studies. 
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Appendix 

A The Superfields 

In addition to the usual MSSM particles the supersymmetric SO(10) contains right-handed 
neutrinos and their scalar partners. With y = x — iOaO we have the following chiral superfields^] 



Hi 


= H i (y) + V2 9H i (y)+99F Hi (y), 


(A.l) 


Q 


= q(y) + V2 6q L (y) + 99F Q (y) , 


(A.2) 


L 


= l(y) + V2 9k(y)+99F l (y), 


(A.3) 


U c 


= Uc(y) + V2eu n c (y) + 99F uc (y) , 


(A.4) 


D c 


= D^y) + V29d K c {y) + 99F DC {y) , 


(A.5) 


E c 


= Ec(y) + V2 9e R c (y) + 99F E c(y) , 


(A.6) 


N c 


= N r c(y) + V29u K c (y) + 99F N c(y). 


(A.7) 



Hi denotes the two Higgs-doublets, 




and 



Ho 



Ht 
Hi 



(A.8) 



Q and L stand for the left-handed quark and lepton doublets and U c , D c , E c and A^ c are the 
right-handed singlet fields. 

Besides the usual bispinors for the quarks and the charged leptons we can introduce 
Majorana-spinors for the right- and left-handed neutrinos 

/ 

and v 



N 



V7° 



(A.9) 



In the symmetric phase of the MSSM no mixing occurs between the fermionic partners of gauge 
and Higgs bosons. Therefore, we have two Dirac higgsinos 



h° 



(MA 



m 



and 



h- 



a 



(A.10) 



2 We are using the conventions of ref. 



with metric 7]^ v = diag(l, — 1, — 1, — 1). 
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which again form an isospin doublet, 



h° 

; (A.n) 

-h~ 

B Boltzmann equations 

The microscopic evolution of particle densities and asymmetries is governed by a network of 
Boltzmann equations. In the following we will compile some basic formulae to introduce our 
notation^. 

It is usually a good approximation to assume Maxwell-Boltzmann statistics, so that the 
equilibrium number density of a particle i is given by 

n?(T) = J d 3 Pl fP with fP {Ek, T) = e~^ T , (B.l) 

where gi is the number of internal degrees of freedom. This particle density can be changed 
by interactions and by the expansion of the universe. Since we are only interested in the effect 
of the interactions it is useful to scale out the expansion. This is done by using the number of 
particles per comoving volume element, 

Yi = ^ , (B.2) 

where s is the entropy density, as independent variable instead of the number density. 

In our case elastic scatterings, which can only change the phase space distributions but not 
the particle densities, occur at a much higher rate than inelastic processes. Therefore, we can 
assume kinetic equilibrium, so that the phase space densities are given by 

m,T) = ^e- E ^ T . (B.3) 

In this framework the Boltzmann equation describing the evolution of a particle number Y$ in 
an isentropically expanding universe reads |14|,9|| 



dY^ z x - 

dz sH (m^) a ^ 



ySyV" ^& + a + ...^i+j + ... 



1p ±a 

V- 1 

B.4) 



YY 

Y e«yc q ' 7 eq (i + 3 + ■ ■ ■ -> 4> + a + ■ ■ ■) 



where z = m^/T and H (nif) is the Hubble parameter at T = m^. The 7 eq are space time 
densities of scatterings for the different processes. For a decay one finds II 



Id ■= 7 eq (^ + ...) = T , (B.5) 



3 For a review and references, see [ H . 
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where Ki and K2 are modified Bessel functions and T is the usual decay width in the rest system 
of the decaying particle. Neglecting a possible CP violation, one finds the same reaction density 
for the inverse decay. 

The reaction density for a two body scattering reads 



7 eq O + a^i + j + ...) 



T 



64vr 4 



(B.6) 



where s is the squared centre of mass energy and the reduced cross section &(s) for the process 
if) + a—*i+j + ... is related to the usual total cross section a(s) by 



o-(s) = - Oty • Pa) ~ m,m a a(s) 



(B.7) 



C Reduced cross sections 

In this section we will collect the reduced cross sections for all the 2 «-> 2 and 2 «-> 3 pro- 
cesses that we had discussed in section [| The corresponding reaction densities, which can be 
calculated analytically in some interesting limiting cases, will be discussed in the next section. 

C.l Lepton number violating processes mediated by a right-handed 
neutrino 

We have mentioned in the main text that we have to subtract the contributions coming from 
on-shell (s)neutrinos, i.e. we have to replace the usual propagators by off-shell propagators 

1 x — a,- ,1 



Dj(x) ' (x-ajY + ajCj 



and 



Dj{x) ' (x-aj) 2 + a jCj 



(C.l) 



To begin with, let us specify the reduced cross sections for the reactions depicted in fig. |3|. For 
the processes I + h I + h and I + H 2 <-> I + H\ one has 



2n I j v ' n x 



x 



+ 



X 



+ 



X 



dj Dj{x) 2D](x 



1 I X + aj 1 In ' X + aj 



Dj(x) 



« j 



X 



+ 



X 



+ 



X 



x [Dj(x) D n (x) Dj(x)D n (x) 



(C.2) 



+ (x + a,- 



/X \ + ^ + Q ^ 



In 



X I ^77, 



D n (x) j v Uj 

where n and j are the flavour indices of the neutrinos in the intermediate state. The interference 
terms with n 7^ j are always very small and can safely be neglected. 
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The reduced cross section for the process I + h <-> / + W 2 reads 



2 dj_ 

33 X 



E (^a, 

3 



_ rf np r f 

tAj Jkj iXj 

+ / T + 



1 



1>3 



x + cij Dj(x) 2D)(x) ' \^ ^j(^) 

2 



In 



x + a 



X 



+ ~7^~~7 — T + 



X 



Dj(x) D n (x) Dj{x)D n {x) 



(C.3) 



— a 



2 1 

+ 



<2 n dj _D n (x) 



111 



X + dj 



2 1 

+ 



DAx) 



In 



37 | ^77, 



The same result is valid for the CP conjugated process. 

~ ~t t 
For the process Z + h <-> I + H 2 one finds 



4 4) (*) 



2^ 



X 



X 2 X 2 X fx + CLj 



+ 



+ 



djlx + aj) D ( x ) Dj-(x) 



+ E Re f( At A 



3<n 



X 



2x 2 



+ X 



\ a 3 

2 1 
+ 



D j (x)D n (x) \a n - aj D n (x) 



In 



X + dj 



+ X 



2 1 

+ 



aj-d n DAx) 



In 



x + a T: 



(C.4) 



rt 



For the scattering I + H 2 — > Z + Z7 C + g and the corresponding CP transformed process we have 



4 5) (*) 



3 a 



57T 



X X 
— + 



+ 



X 



+ E Re [(At A,) 



^JdrJ- 



X X 

+ =^7T + 



" 1 + 



x + a,- \ / x + a 



x 



^•(x) 



In 



a ; 



x 



(C.5) 



+ (x + a.,) 



v On-% £> n (x) 



111 



x + a,- 



a ; 



+ (x + a n ) 



In 



Finally, we have two processes which do not violate lepton number but merely transform leptons 
into scalar leptons and vice versa. We have the 2 — > 2 scattering Z + H 2 «-> Z + h, 



^ 6) (*) = ^E (At A,) 



47T 4- 



"J 



2 x 2 



Dj(x)D n (x) 



(C.6) 



and the 2^3 process Z + Zt <-> Z + + [/ 

3a u 



167r2 In 



E (AtA„) 



2 x 2 



Dj(x)D n (x) 



Let us now come to the 2^3 processes shown in fig. |j. 
For the transition q + U c ^l + l + H 2 the reduced cross section reads 



(C.7) 
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4 8) ( 
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3a u 
16^2 



~ + ~T=^ 

a 3+ c j Ja jCj 



arctan 



CljCj 



+ arctan 



In 
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The remaining integral cannot be solved exactly. However, it can be neglected for x > dj,d n 
and for x < dj, d n it can be approximated by 



1 



dxi r= ln 

CL j | $"t% ^ 



dj + djCj 



(C.9) 



ln 



Inf^Uspf * Uspf an ~ X ) , 

\ d n J \d n + dj-x) \d n + dj-xj 



where Sp(x) is the Spence function. 
For the scatterings q + U 
cross sections are equal, 
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~t ~ ~ ~t 

For the process I + q —> I + U c + H 2 and similar reactions one gets 
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The remaining integral can again not be solved exactly. However it can be approximated by 
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for x < d n and for x > d n it can be neglected. 

Finally, we have to compute the t- and u-channel processes from fig. [5] which give compara- 
tively simple contributions. 

For the processes I + I <-> h + h and I + I <-> H\ + H\ we get 
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At this order of perturbation theory the same result is valid for the CP transformed processes. 
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For the scattering I + I <-> h + H\ one has 
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The 2-^3 process H 2 + q ] ^ l ] + l ] + U c gives 
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For the related transition I + I ^ U c + q + H 2 we have 
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There are some 2^2 processes left which do not violate lepton number but simply transform 
leptons into scalar leptons, like in the process 1+ I *-+h + H 2 
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or in the similar process I + H\ <-> h + I 
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Finally, the last process i + r ^ h + q^ + u c gives 
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C.2 Scattering off a top or a stop 

For the processes specified in fig. || the reduced cross sections read 
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To regularize an infrared divergence in the t-channel diagrams we had to introduce a Higgs-mass 
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In the calculations we have used the value \x = 800 GeV. 

The analogous processes involving a scalar neutrino (cf. fig. 0) give similar contributions 
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C.3 Neutrino pair creation and annihilation 

With the abbreviations 
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the reduced cross sections for the right-handed neutrino pair creation read 
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For the scalar neutrinos one has similarly 
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For the diagrams involving one neutrino and one sneutrino (cf. fig. [TO]) one finally has 
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D Reaction densities 



In general the reaction densities corresponding to the reduced cross sections discussed in ap- 
pendix |C| have to be calculated numerically. However, there exist some interesting limiting 
cases where one can calculate them analytically. 



D.l Lepton number violating scatterings 

In the Boltzmann equations we do not need every reaction density 7 
(cf. sect. |3.1| ). We only have to consider the combined reaction densities 
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For low temperatures, i.e. z ^> lfy/a] , the dominant contribution to the integrand of the 
reaction densities comes from small centre of mass energies, i.e. x <C dj. In this limit the 
reduced cross sections a N for the L + L violating or conserving processes behave differently. 

(D.5) 



For the L + L violating scatterings (i — 1, . . . , 5, 8, 10, 12, ... , 16) one finds 

af (xx for x <C aj , 
while one has 
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for the L + L conserving processes (i = 6, 7, 9, 11, 17, 18, 19). Hence, the reaction densities can 
be calculated analytically in this limit and one finds 
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For high temperatures, i.e. z <C l/^/oJ , we can use the asymptotic expansions of the reduced 
cross sections to compute the reactions densities and we get 
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where 7 E = 0.577216 is Euler's constant. These reaction densities are therefore proportional to 
T 4 at high temperatures, as expected on purely dimensional grounds. 

For intermediate temperatures z ~ the reaction densities have to be computed 

numerically. This becomes increasingly difficult in the narrow width limit, where 1/Dj(x) 
has two very sharp peaks. However, in the limit Cj — > the two peaks in 1/Dj(x) cancel each 
other, since they have a different sign, while the peaks in 1/ D 2 -(x) add up. Therefore, the terms 
proportional to 1/Dj(x) or l/Dj(x)D n (x) with n ^ j can be neglected in the narrow width 
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limit, while l/D 2 Ax) can be approximated by a 5-function 
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An analogous relation holds for 1/Dj (x). 

These relations allow to calculate the contributions from the s-channel diagrams to the 
reaction densities analytically in the limit Cj — > 0, while the contributions from the t-channel 
diagrams can easily be calculated numerically. 

D.2 Interactions with quarks and squarks 

The reaction densities 7^ for the interaction of a (s)neutrino with a top or a stop can also 
be calculated analytically in the limit of high temperatures z <^ 1/ ^/a] . For the s-channel 
processes one finds 
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For the t-channel reaction densities one has analogously 
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7^ and 7^ are several orders of magnitude smaller than the other 7^ for small z and can 
therefore be neglected at high temperatures. 

By using the series expansions of the Bessel functions, one sees that the processes with a hig- 
gsino in the t-channel, i.e. 7^, 7^, 7^ and 7^, behave like T 4 ln(T/Mj) at high temperatures, 
whereas the other reaction densities are proportional to T 4 . 

D.3 Pair creation and annihilation of neutrinos 



In the Boltzmann equations we only need certain combinations of reaction densities which can 
easily be evaluated for high temperatures: 
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i.e. these reaction densities are proportional to T 4 ln(T/(Mj + Af,-)) at high temperatures. 
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